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We calculate the neutrino emissivity of superfluid neutron matter in the inner crust of neutron 
stars. We find that neutrino emission due to fluctuations resulting from the formation of Cooper 
pairs at finite temperature is highly suppressed in non-relativistic systems. This suppression of the 
pair breaking emissivity in a simplified model of neutron matter with interactions that conserve 
spin is of the order of Vp for density fluctuations and Vp for spin fluctuations, where vp is the 
Fermi velocity of neutrons. The larger suppression of density fluctuations arises because the dipole 
moment of the density distribution of a single component system does not vary in time. For this 
reason, we find that the axial current response (spin fiuctuations) dominates. In more realistic 
models of neutron matter which include tensor interactions where the neutron spin is not conserved, 
neutrino radiation from bremsstrahlung reactions occurs at order Vp. Consequently, even with the 
suppression factors due to superfiuidity, this rate dominates near Tc- Present calculations of the 
pair-breaking emissivity are incomplete because they neglect the tensor component of the nucleon- 
nucleon interaction. 



I. INTRODUCTION 



The long term (10^ — 10^ yrs) cooling of isolated neutron stars and the thermal evolution of accreting neutron stars 
in binary systems are sensitive to the weak interaction rates at high density. A recent review of the theory, modeling 
and observational constraints on neutron star thermal evolution can be found in Ref. [ij. Here, we focus on a specific 
neutrino process which is expected to be relevant in superfluid neutron matter. 

In superfluids, Cooper pairs break and recombine constantly at finite temperature. Such processes can dominate the 
density and spin-density fluctuations. In pioneering work, several decades earlier. Flowers, Ruderman and Sutherland 
recognized that these fluctuations can couple to neutrinos through the weak neutral current. They showed that when 
the temperature was less than but comparable to the critical temperature for superfluidity, neutrino emission due 
to the Cooper pair recombination processes was important 0. Subsequently, this process which is now commonly 
referred to as the pair-breaking and formation (PBF) process and was recomputed in Refs. [1, 01 ■ Its role in the 
thermal evolution of isolated neutron stars was studied and shown to important 0, 0, 0] ■ 

Another context in which the PBF process plays a role is in accreting neutron stars which exhibit x-ray bursts and 
superbursts. Current models for superbursts indicate that they arise due to unstable burning of carbon in the ocean 
of accreting neutron stars Q. Agreement between theoretical models (for the light-curves and recurrence times) and 
observation relies on the assumption that carbon is ignited at a column-depth of about 10^^ g/cm^. However, to ignite 
carbon at this depth the temperature there should be ~ 5 x 10* K [1, [lOl- Consequently the ignition condition and 
hence the recurrence time for superbursts are sensitive to the temperature profile of the crust in accreting systems. 
The temperature profile in turn depends on the balance between heating in the crust due to electron captures and 
pycnonuclear reactions [ll| and neutrino cooling Q . 

The inner crust of the neutron star is expected to contain a neutron superfiuid. When the matter density exceeds 
the neutron drip density (4 x 10^^ g/cm^) a relatively low-density neutron liquid coexists with a lattice of nuclei p^ . 
Here, the attractive s-wave interaction between neutrons induces superfluidity. The pairing-gap rises from zero at 
neutron drip to a maximum of about 1 MeV when the neutron fermi momentum /cf ~ 200 MeV and then decreases 
to zero in the vicinity of the crust-core interface. For temperatures of relevance to the accreting neutron stars, it 
was found that the PBF process in the neutron superfluid resulted in rapid neutrino losses and cooled the crust to 
temperatures below those required for carbon ignition at the favored depth p^ . Subsequently, additional heating 
processes in the outer crust due to electron captures on nuclei were shown to be relevant but were unable to produce 
the necessary heating in models which included PBF process in the crust [l^ . 

The preceding discussion motivates a detailed investigation of the neutrino emissivity arising due to the PBF process 
in neutron matter in the inner crust. Recently, this was recalculated by Leinson and Perez who find that earlier 
calculations violated vector-current conservation p^ . [l5j . An improved treatment that satisfies current conservation 
yielded a result that was suppressed by the factor Vp where vp = kp/M is the neutron fermi velocity, kp is the neutron 
fermi momentum and M is neutron mass. In the neutron star crust, where vp 0.1 this suppression is significant. 
Subsequently, Sedrakian, Muther and Schuck also calculated the PBF rate using an improved treatment based on 
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Landau Fermi liquid theory and found that it was suppressed by the factor ~ T/M where T is the temperature 
p^ . This suppression is parametricaUy different from that obtained in Ref. The PBF rate was also recently 

examined in Ref. [l3| where the authors also accounted for Fermi liquid effects in superfluid neutron matter using the 
Larkin-Migdal-Leggett formalism [3, [3 ■ They find that the vector current response is suppressed by the factor v% 
in agreement with the finding of Leinson and Perez. 

In this article, we reexamine the nature of density and spin-density fluctuations in superfluid neutron matter using 
a simplified nuclear Hamiltonian. Our main findings are 

1. The spectrum of density fluctuations is suppressed at order in pure neutron matter in agreement with the 
findings of Leinson and Perez. 

2. The u|i suppression is not generic and is not a consequence of vector current conservation. It is specific to simple 
one component systems where all particles have the same (weak) charge to mass ratio. In multicomponent 
systems such as the neutron star crust where neutrons coexist and interact with nuclei, the density fluctuations 
of the neutron superfluid occur at order . 

3. For the case of simple nuclear Hamiltonian with only central interactions that conserve spin, spin-density 
fluctuations occur at order and the these fluctuation dominate the neutrino emissivity. 

4. For the case of realistic nuclear interactions which contain a strong tensor component, spin is not conserved, 
and spin fluctuations arise at order w^. This feature is well known in the context of neutrino emission from 
neutron-neutron bremsstrahlung. We find that the bremsstrahlung rate continues to be the dominant neutrino 
emission mechanism even in the superfluid state for T > Tc/5. 

The article is organized as follows. We begin by discussing the relation between the neutrino emissivity and the 
density and spin-density response functions. This is followed by a detailed investigation of the density-density response 
function and the role of vertex corrections in the superfluid state. Here we show that the vertex corrections required 
by conservation laws strongly suppress the response relative to the predictions of mean field theory as suggested in 
earlier work. Subsequently we discuss the spin-density response function and show that it dominates over the density 
response. Finally, we will discuss the various contributions to the neutrino emissivity and conclude that neutron- 
neutron bremsstrahlung rate is typically larger PBF process even in the superfluid phase. We conclude with a critical 
discussion of our study here and related earlier work. We recognize that all calculations of the neutrino rates from 
PBF are missing a key aspect of the nuclear force - namely the tensor interaction. 



II. NEUTRINO EMISSIVITY AND RESPONSE FUNCTIONS 



The neutrino emissivity is defined as the rate of energy loss per unit volume and is given by 

ivv - ^ 



where Gp is the Fermi weak coupling constant, fc = {uj,k), gi=i,2 are the on-mass-shell four-momenta of neutrinos, 
L"^{qi,q2) = Tr [7^(1 — 7^) ^i7''(l — 7^) ^2] and H^f^iq) is the retarded polarization tensor [l^. Using Lenard's 
identity [2^ , we can simplify Eq. [1] to obtain 

fy2 r poo 

192 TT^ 7 Jo exp (/3w) - 1 

where the superfluid response function i?Q^(a;,|fc|) in general contains both the vector and axial- vector response 
functions and is given by 

R^p{~u:, \k\) = -cl 5m[nj:^(c., |fc|)] - c\ 5m[n^^(a;, |fc|)] . (3) 

In the non-relativistic limit, we shall focus on density fluctuations and ignore velocity fluctuations. In this case the 
vector-polarization function H^^(a', |fc|) ~ 5'^ 5^^ IIo(a;, |fc|) where Ho(w, k) is the density-density correlation function 
(2H given by 

Ho(c., |fc|) = -i [ d^x e-^^^-^-^'^^TripG \p{x, t), p(0, 0)]) , (4) 
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where pa is the density matrix and p{x, t) is the density operator. Similarly the axial response in the non-relativistic 
limit is dominated by spin fluctuations and we can write n^^(ci;, |fc|) = (5^ (5^ ny (ci;,|fc|) where i,j = 1,2,3 and 
njj(w, k) is the spin correlation function [2l| given by 



n„ (a;, |fc|) = -i / d^x e-'(^-^-^'^Tr{pG t), a,(0, 0)]) . (5) 



The diagonal components in Eq. Oare equal and is denoted by Her, while the off-diagonal components of Ily do not 
contribute to the emissivity of an isotropic medium We can therefore write the neutrino cmissivity as 



192 



jd?k e [cl Ip{k) + ic\ I^{k)\ , (6) 



where 



Ip{k) = - do; ^^<^m[Uoiu,\k\)], (7) 

Jk exp {f3uj) - 1 

lAk) = - Tda; f^^-|) 3m[n.(c.,|fc|)]. (8) 

Jk exp(/3w)-l \k^ 3 J 

III. VECTOR RESPONSE 

First, we calculate the vector-current response function to verify and understand the nature of suppression factors 
found in Refs. [13, [3 ■ Calculations of the superfluid density response has a long history in condensed matter physics 
and a pedagogic discussion can be found in Ref . [2^ . Here we describe neutron matter at low density with the model 
Hamiltonian 

'^^Yl ^pOfcT^fcT +^Xl"feT"-fci"'^Ta-fci, (9) 

p,spin p,p' 

where V is the effective four-fermion interaction. We regulate the short-range interaction by using a momentum 
cut-off. The re- normalization scheme is implemented by specifying the gap and requiring that V{A) satisfy the gap 
equation at zero temperature 







(27r)3 2Ep 



To describe neutron matter at densities of relevance to the crust, we choose to display results at kp/{2m) = /.t = 30 
MeV and a momentum cut-off of A = 2kp. Though we are ultimately interested in calculating the PBF rate in the 
temperature range T ~ 1 — 10 x 10* K and near kp < 100 MeV where the gap A < 0.1 MeV [2^, our numerical 
results are primarily at slightly larger densities and for A = 1 MeV where the numerical computations are somewhat 
easier. Although finite-range effects of the nuclcon-nuclcon interaction are relevant in computing the magnitude of 
the pairing gap in the neutron star crust here we will restrict our analysis to a simple zero-range interaction but a 
strength adjusted to reproduce the pairing gap in more sophisticated calculations [2J|. 
We define the polarization tensor in the mean field approximation by 

nr/3(^, 1^1) = / (^Tr [7„G(p + fc)7^G(p)] , (11) 

where 7^ = (^t^, 1 (p + k/2)/M^. Here the zero-zero component of H^^ corresponds to the density-density response 
function defined in Eq. [D Explicitly this is given by 

Hmf(c., \k\) = -I j ^Tr [t^G{p + k)T^G{p)] , (12) 

where the quasi-particle propagator in the superfluid state is given by 



(13) 
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and 1 is the unit matrix, Ti=i.2,3 are the 2x2 Pauh matrices (acting in the Nambu-Gorkov space). The quasi-particle 

energy is Ep = + IS? where = {p^ /2m — fi) and A is the superfluid gap [2^ . 

The mean-field polarization tensor violates current conservation and the F-sum rule. This is a well established 
finding and a lucid discussion can be found in the original papers by Anderson (25j and Nambu (26| . To restore 
current-conservation it is necessary to replace the bare vertex function 7q(p -I- q,p) by the dressed vertex T^(p + 
k,p) corresponding the dressed quasi-particles in the superfluid. The dressed vertex is then required to satisfy the 
generalized Ward-identity (GWI) given by 

a.r"(p + k,p)- k,r{p + k,p)= G-\p + k)T3 - T3G~\p) , (14) 

where j = 1, 2, 3 and as before the four-vector k = {lo, k). However the GWI does not uniquely determine the vertex 
function. It must be obtained explicitly as a solution to an integral equation which describes the modification of the 
vertex from the medium. This integral equation for the vertex, obtained in the random-phase approximation (RPA), 
is known to satisfy the GWI and is diagrammatically in Fig[T] Explicitly this is given by 



Ta = 7a + i ^ / jtAi rs G{q + k) r„ G{q) T3 . (15) 
J (27r)^ 

We are interested in the zeroth-component of Eq[Tn] since it is this that affect the density response. In this case an 




FIG. 1: Diagrammatic representation of the RPA vertex equation. The wavy line is the external weak current and curly line 
represents the short-range strong interaction which in this work is simplified to a point interaction. Solid lines represent the 
quasi-particles and the dark circle is the dressed vertex. 



approximate solution has the form 



r 1 2k . 

r, = ^r, + ^^r, (16) 



In weak coupling, x — ^ -^(0) 1 where V is the four-fermion coupling and A^(0) oc A/ kp /tt \s the density of states 
at the Fermi surface and can be neglected. However k has an essential singularity at T = corresponding to existence 
of a Goldstone excitation that couples to density fiuctuations. At T = and weak coupling 

(17) 



where speed of the Goldstone mode Cs — kp/^/iM . This approximate form for the vertex was used in Ref. [lB| to 
compute the response function. It is possible to solve the vertex equation (Eq. [T5)) to obtain the exact solution in 
the random phase approximation [27| . In Fig. [5] we compare the exact solution (solid- line) to the vertex equation 
to the approximate expression above (dashed- line). For a better approximation to the full vertex, we can employ a 
shifted value of by shifting it so that the approximate and full vertices diverge at exactly the same value of go- 
This "shifted" vertex is also plotted in Fig. [2] with a dotted line, but the result matches sufficiently well with the 
full vertex that it is not visible. We also find that this relative agreement between the real parts of the full and 
approximate results is not strongly modified at finite temperature. In the region where lo > 2A where we are above 
the threshold for producing quasi-particles we can expect that the vertex equation will have a non-zero imaginary 
part. The approximate solution to the vertex equation in Eg. 1171 neglects this contribution. The imaginary part of 
the full vertex is plotted in the right panel of Fig. [21 As expected below ui = 2A, the imaginary part of the response 
vanishes as required from the delta function given by the imaginary part of the propagator in Eq. 1181 The magnitude 
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FIG. 2: Real parts of the dressed vertex function at fixed momentum transfer. The approximate solution in Eq.[T7]is in excellent 
agreement with the "exact" RPA vertex especially when the speed of sound is shifted to match the pole structure. Right panel 
shows the imaginary part of the vertex function is finite but small. The approximate vertex assumes that imaginary part is 
zero. 



of the imaginary part is much smaller than the real part, but as we shall discuss later could make an be important 
contribution to the response. 

To help make contact with earlier results obtained in Ref. 0, [3 we discuss the response function in different 
approximations. First, wc obtain the mean field response at T = by doing the po integration in Eg. 1121 



1 1 



Iq 



lo = ^ — -IT^ (1^) 



In the long-wavelength limit {k 0) this can be simplified further and wc find that 



Umf{u;, \k\ ^ 0) = y ^ lo . (19) 

For LO > 2A and fc = the mean field result predicts a non-vanishing response given by 

5m[nMF(a.)] = f ^4^_ + ^) . (20) 

As mentioned earlier a non-zero response at /c = violates current conservation and the related F-sum rule given by 

cLu uj 3m[nMF(c^)] = {[[H,p{k)],p{k)]) , (21) 

where the RHS vanishes in the k = limit when p commutes with the Hamiltonian. This is simply a consequence of 
the well known fact that radiation can arise only as a result of particle acceleration for conserved charges. The vertex 
correction discussed earlier remedies this problem. The RPA response function is defined by 

nRPA(t^, \k\)^-ij Tr [r^Gip + k)T^G{p)] , (22) 

where To is the dressed vertex that satisfies Eq. 1151 First we use the approximate form of the vertex given in Eq. [17] 
to obtain the zcro-tempcrature density-density polarization function. As wc showed earlier the approximate vertex 

2Au . 

To = 7-3 -5 37:5 IT2 (23) 
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provides a very good description of the real part of the exact result. Substituting this into Eg. 1221 and performing the 
integration over the po we obtain 

^^^^ ' ' J 2(2^)3 1^ EpEp+, J " ^ ^ 

The response function is related to the imaginary part which is explicitly given by 

Qm[n^M^,k)]^-J^JSico~E,~E,+,) where J = (^1 - M£±^^^±^) (25) 

It is straightforward to verify that the imaginary part of IIrpa vanishes in the limit k —f as required by current 
conservation. At finite temperature and for the full vertex this continues to hold and the response vanishes at fc = 
because of the gap equation. The critical question then is to inquire if the order fc^ term also vanishes. To address 
this we expand the zero-tempcrature RPA response in Eq. [25]in powers of k. Expanding the integrand in Eg. [25l we 
obtain 



J = 



2E^ 
p 



? 2 



+ O [c, kf (26) 



where x is the cosine of the angle between the momenta p and q. In the long-wavclcngth most of the support from 
the energy delta function is near p = pp. Further, in this region we can approximately replace x^ by its mean value 
given x^ ~< x^ >= 1/3. This implies that the quadratic term nearly vanishes since Cs = pp/i'm^/S). The nature of 
this cancellation depends on the value of Cg employed in the approximated vertex. When we use the self-consistent 
RPA vertex obtained by solving Eq. [12] we indeed find that this cancellation is nearly exact and Eq. [H] receives 
contributions at order fc^. 

The imaginary part of the response function at zero temperature obtained in different approximate schemes are 
shown in Fig. [3l The striking feature is that results obtained using the dressed vertex are suppressed because 
they are oc and the response using the self-consistent RPA vertex labeled "Full RPA" in the figure is even more 
suppressed because the relevant contribution occurs at order cf. The result labeled "Approx." is obtained using 
Cs = Cs,o = Pf/ {V^M). We note that the density response function calculated within RPA coincides with earlier 
calculations by Kundu and Reddy Ref . [131 [Ml • From the Fig. [3] it is clear that the magnitude of the suppression 
obtained using the approximate vertex is sensitive to Cs- This is because of the large cancellation at order k^. To 
understand the nature of this cancellation we first note that vector current conservation docs not require the order 
Vp fc^ contribution to vanish. This would require fine tuning the velocity of the Goldstone mode to precisely cancel 
the contribution at order fc^. In general, we are aware of no mechanism that can accomplish this in a multicomponcnt 
system. 

The preceding arguments raises the following question: what is the underlying physics responsible for the large 
cancellation at order k^ occurring in the RPA calculation of the response function for pure neutron matter ?. To 
address this question we note that a similar cancellation at order fc^ occurs in neutrino emission due to density 
fluctuations in the bremsstrahlung reaction nn — > nnvv in the normal phase. Here the square of the matrix element 
contributes only at order v%k'^ [H, [2^. This result is well understood, especially in the context of electromagnetic 
bremsstrahlung in proton-proton collisions where radiation occurs due to time variation of the quadrupole moment. 
The dipole radiation vanishes because the dipole moment for identical particles interacting with each other does not 
vary in time. In fact for any system of N identical particles with charge e and mass m, the dipole moment 

g g 

d = S^efi = —T,im fl = —Y^im Rqu (27) 
m m 

does not vary because momentum conservation ensures that Rcm = 0. In the quantum mechanical calculation of the 
bremsstrahlung process this arises due to the destructive interference between the amplitudes for radiation from the 
two charges. This interference ensures that the square of the matrix element vanishes at order vl^ k^ and the leading 
contribution is quadrupolar and occurs at order v^^ fc^. 

When the (weak) charge to mass ratio of the two particles is different, the cancellation at order w|, k^ would 
be absent. For example bremsstrahlung from neutron-proton scattering occurs at order v'^ k^. Hence we argue in 
multicomponcnt systems there is no symmetry that requires the sensitive cancellation at this order. In the neutron 
star context, where neutrons in the crust couple to a background lattice of neutron-rich ions, we can expect that these 
interactions can induce a response at order k^. In our calculation, such interactions can for example induce a shift 
in the superfluid velocity Cs by the polarization of the lattice. As we have seen in Fig. [3l even a small shift leads to 
a relevant contribution at order c^fc^. Thus we conclude that in realistic situations a non- vanishing density response 
at order cj.k^ is to expected. However, as we show below this response is still small compared to the axial-current 
response. 
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FIG. 3: Imaginary part of the density-density polarization at T = obtained using different approximations. The vertex 
corrections induce a large suppression in the response relative to the mean-field result. The result obtained using the self- 
consistent RPA vertex (shown with the dotted line) is suppressed by a factor ~ 10~*. See text for a discussion relating to the 
response obtained using the approximate vertex. 



IV. AXIAL CURRENT RESPONSE 



111 the non-relativistic limit, the diagonal part of the axial polarization tensor defined in Eq. [S]can be written in 
terms of the Nambu-Gorkov propagators [23|. To begin we will ignore vertex corrections because there is no Goldstone 
mode associated with spin-fluctuations in the case of singlet pairing and discuss the one-loop mean field polarization 
tensor. In this case we can write 



n,(c^,|fc|) = - 

At T = we can do the pq integration to obtain 



Tr [i Gip + k) i Gip)] 



d^p 
2(27r)3 

and the imaginary part in the region > is given by 



Qm[U^{uj,\k\)] = - 



d^p 
2(27r)3 



(28) 



(29) 



(30) 



It is easily verified that the imaginary part of Eg. 1301 vanishes at fc = 0. Expanding the integrand in Eq. [2D]in powers 
of k we find that 



9m[n,(c^,|fc|)] 



1 



327r2 



dpp^ I dx [J^^x" e + 0[k^\ ) 5{u: -Ep~ Ep+k) 



m? E^ 



(31) 
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where x is the angle between k and p. In the long-wavelength limit the delta function provides support only in region 
p Pp. Here we can further simply the result to obtain 



Sm[no-(ti', \k\ 



1 



487r 



2 7 2 

Vp k 



2 ""F 



dp P ^ dxS{uj-Ep-Ep 



+ k) 



0[k]^ 



(32) 



In Fig. [4] we plot the axial response function. The vector response obtained in different approximations discussed 
earlier is also shown for comparison. The results indicate the axial response is significantly larger because of the large 
cancellation in the vector response at order k^. In situations where this cancellation does not occur, as in the case 
when we set ~ O.SCg^Oj the vector and axial response functions can become comparable. 
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FIG. 4: The imaginary part of the axial response function in the non-relativistic limit. 

The one- loop axial response in Eg. 1281 vanishes in the k ^ limit. This is consistent with the F-sum rule for the spin 
response associated the model Hamiltonian in Eq. [9l However it is well-known that the realistic nuclear Hamiltonian 
which contains both tensor and spin-orbit interactions does not commute with the spin operator. This implies that 
for realistic nuclear interactions which include pion exchange 

/oo 
duj Lu\tk-.o^m[n,{uj)]^ {[[H,a{k^O)la{k = Q)]) ^0. (33) 
-oo 

Consequently, the correct axial response obtained using realistic interactions has to be finite in the limit fc = 0. 
From Eq. [32]we see that the mean field response function violates this expectation. Corrections to the one- loop axial 
response function arising due to non-central interactions such as pion exchange are therefore will therefore be critical 
in the long- wavelength limit. 

The importance of pion exchange in the axial response and neutrino emissivity was already realized in the pioneering 
work of Friman and Maxwell [28| . Although their calculations were only applicable to the normal phase of neutron 
matter, they showed that pion exchange was more important than central nuclear interactions. In particular, they 
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found that the matrix element for neutrino emission from neutron-neutron bremsstrahlung in the axial channel was 
non-zero in the long-wavelength limit when they included one-pion exchange. It is therefore critical to incorporate 
these non-central interactions in the calculation of the axial response function of the neutron superfluid. This would 
require that we include both vertex corrections and two-loop effects which include 2p — 2h excitations in the superfluid. 
This work is in progress and will be reported elsewhere. 



V. NEUTRINO EMISSIVITY 



The neutrino emissivity for the processes discussed above is plotted in FigO including the axial and vector contri- 
butions to the PBF emissivity and the bremsstrahlung emissivity with and without the suppression factors resulting 
from superfluidity as obtained in Ref. [3| ■ The emissivity from the axial response dominates and is two or three order 
of magnitudes larger than the full RPA vector response. A shift of the speed in sound of the Cooper pair excita- 
tions will modify the emissivity quadratically in Cs as expected from the imaginary part of the vector response. The 
bremsstrahlung emissivity is larger than the pair-breaking emissivity except at low temperatures when the suppression 
of the bremsstrahlung from superfluidity is strong. 
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FIG. 5: The neutrino bremsstrahlung and pair-breaking emissivities. The dashed line is the bremsstrahlung emissivity in 
non-superfiuid neutron matter and the short dashed-dotted hne is the corresponding emissivity in superfluid neutron matter. 
The long dashed-dotted line is the axial part of the pair-breaking emissivity and the dotted line is the vector part of the 
pair-breaking emissivity using the full RPA vertex. The soHd line is the vector pair-breaking emissivity using the approximate 
vertex. 



VI. CONCLUSIONS: 



We have studied the one-loop vector and axial current response functions of relevance to neutrino emission in 
superfluid neutron matter. Through explicit calculations we found that there is strong suppression of the vector 
resp onse when vertex corrections are included. This is in close agreement with the findings of Leinson and Perez 
For pure neutron matter, the RPA vertex function does indeed show that the supression factor is of order Vp. 
However we have shown this suppression arises not only because of vector current conservation, but also because the 
radiation in the vector channel for pure neutron matter occurs only due to time variation of the quadrupolc moment. 
In the realistic context where neutrons interact with the background lattice of of neutron-rich ions the suppression 
in the vector channel is of order Vp. We showed that even a small shift in the speed of the Goldstone mode due 
to the lattice can make a relevant contribution to the response at this order. Finally, we showed that both of these 
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emissivities are likely smaller than the neutrino bremsstrahlung emissivity (including the suppression factors from 
superfluidity) unless the temperature is significantly smaller than the critical temperature. 

The axial response function was shown to be numerically more important because because the dominant contribution 
occurs at order Vp. Although this was the relevant contribution in our model, we demonstrated that tensor interactions 
arising due to pion exchange would lead to important corrections to this estimate. The F-sum rule in this case indicates 
that these correction would result in a non- vanishing response at order Vp. In the normal phase the tensor interaction 
is responsible for the emissivity at long-wavelength and we suspect that this will continue to be the case in supcrfluid 
matter. We anticipate that tensor interactions will also affect the PBF fluctuations in the spin channel and can modify 
our results and the regime in temperature where PBF can dominates over the bremsstrahlung emissivity. This issue 
is currently being investigated and will be reported elsewhere. Although it is now clear that the neutrino emissivity 
due to density fluctuations arising from PBF processes in the vicinity of the critical temperature is not significant, our 
present understanding of neutrino processes in the supcrfluid phase remains rather incomplete and warrants further 
study. 
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